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Higher order QED corrections to deep inelastic scattering 
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We calculate the leptonic 0(a 2 L) QED corrections for unpolarized deeply inelastic ep scattering using mixed 
variables. 


1. INTRODUCTION 

Deep inelastic electron-nucleon scattering allows 
for fundamental QCD tests investigating the scal¬ 
ing violations of structure functions in the pertur¬ 
bative regime of large values of Q * 2 . The detailed 
knowledge of the structure functions enables to 
study various aspects of the dynamics of non- 
Abelian gauge theory, and is necessary for the 
future experimental search for the Higgs-boson 
and new particles at TEVATRON and LHC. 

One of the major goals of the experiments HI 
and ZEUS at the ep-collider HERA at DESY is 
to perform a QCD test at large space-like virtual¬ 
ities Q 2 at high precision. This presumes to know 
the QED radiative corrections to the double¬ 
differential scattering cross sections of deeply in¬ 
elastic ep corrections as precisely as possible. Pre¬ 
vious calculations of the radiative corrections for 
the unpolarized cross sections at leading order [1- 
5]@, the leading-log level [7 12] to leading and 
higher orders, and QED-resummations of small¬ 
er terms [11,13] revealed that these corrections are 
very large, in a wide kinematic range of x and 
Q 2 . The corrections are, moreover, complicated 
by a new type of sizeable contributions as the 
Compton-peak [ [hi. This makes it necessary to 
extend the calculations to higher orders. 

The higher order leading-logarithmic contribu¬ 
tions O [(aL) fe ] to QED corrections are obtained 
as the leading order solution of the associated 
renormalization group equations [ 0 for mass 
factorization. These corrections are universal, 
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2 For QED corrections to polarized lepton scattering off 
polarized nucleons see [ kJ . 


process-independent w.r.t. their structure, and 
are given in terms of Mellin-convolutions of lead¬ 
ing order QED splitting functions. The next-to- 
leading order (NLO) corrections can be obtained 
along the same line. However, besides the split¬ 
ting functions to NLO also the respective process- 
dependent Wilson coefficients and operator ma¬ 
trix elements in the on-mass-shell (OMS) scheme 
contribute. In the past this method was applied 
to calculate the 0{a 2 ) initial-state QED radia¬ 
tive corrections to e + e _ —■> p + p~ in Ref. [ |To| . 

In the present paper we summarize results of a 
recent calculation [0 of the leptonic QED cor¬ 
rections to deeply inelastic ep scattering defining 
the double-differential scattering cross section for 
mixed variables [8,2,5] to 0(a 2 L). If compared 
to the case dealt with in [ [hJ to 0(a 2 L) the 
present calculation is more complicated due to 
the emergence of final state radiation, the double¬ 
differential cross section and the relevant rescal¬ 
ing, which implies non-Mellin type convolutions 
in general. We first summarize main kinematic 
aspects and present then the different contribu¬ 
tions to the leptonic NLO QED corrections. 

2. Mixed variables 

ep collider experiments allow to measure the 
kinematic variables defining the inclusive deep- 
inelastic scattering cross sections in various ways 
since in principle four kinematic variables are 
available with the energies and angles of both 
the outgoing lepton and the struck-quark, see 
e.g. [ ||^]. At the Born level all methods are 
equivalent, however, resolution effects as a conse¬ 
quence of the detector’s structure, differ in certain 
kinematic regions. The Bremsstrahlung-effects 



of QED radiative corrections change this picture 
drastically and the QED correction factors de¬ 
pend on the way the kinematic variables, as e.a. 
Bjorken-p and the virtuality Q 2 are measured. u 
In the present paper we calculate the NLO ra¬ 
diative corrections in the case of neutral current 
deep-inelastic scattering for mixed variables, i.e. 
that Q 2 = Q 2 is measured at the leptonic and 
y = yh is measured at the hadronic vertex, and 
x m = Q 2 /(Sy h ). The Born cross section for 7 - 
exchange is given by : 


riV°> 

dydQ 2 
with 


2'kcx 2 r 0 , . „ 

^rhV2xF, + 2(1-y)F 2 ], 


(1) 


F 1 (x,Q 2 ) = ^2 [qk{x,Q 2 ) + q k (x,Q 2 )\ , (2) 

Z k= 1 

F 2 (x,Q 2 ) = 2xF 1 (x,Q 2 )+F L (x,Q 2 ) . (3) 


Here, Fi^,l(x,Q 2 ) denote the nucleon structure 
functions for photon exchange, and q(x,Q 2 ) and 
q(x,Q 2 ) are the quark and antiquark distribu¬ 
tion functions. The sub-system variables obey 
the following rescaling relations for initial- and 
final-state radiation : 

ISR : y=— Q 2 = zQ 2 , S = zS, x = zx m , 
z 

J\z) = 1 , 4 = min {vh, j , (4) 

FSR : y = y h , Q 2 = ^, S = S, £=—, 
z z 

J F (z) = —, Zo=x m . (5) 

z 

Here, J I,F (z) are the initial- and final-state Ja- 
cobians d 2 (y,Q 2 )/d 2 (yh,Qf), and z 0 marks the 
lower bound of the sub-system rescaling variable 
2 e [z 0 ,l]. The rescaling in Eqs. ([|,|]) was cho¬ 
sen such, that both the initial- and final-state 
operator matrix elements can be expressed with 
a variable 2 e [0,1]. Q g is introduced as a scale 
to cut away contributions of the Compton peak. 
Although these terms do formally belong to the 
QED radiative corrections, they stem from a kine¬ 
matic domain of low virtualities and are therefore 

3 See e.g. Refs. [8,2,5] for a comparison a wide range of 
different choices of measurement. 


not being associated to deep inelastic scattering. 
The scale Q q can be chosen by experiment ac¬ 
cepting only those events in the sample to be an¬ 
alyzed for which the hadronic Q 2 is larger than 
Qq. By this measure the Compton peak is cut 
away widely and the QED-correction factor is 
dominated by a deep-inelastic sub-process by far. 
In the case of mixed variables the leptonic QED 
radiative corrections can be easily grouped into 
those for the initial and final state. The separa¬ 
tion scale between the two kinematic regions is 

Ql 


3. NLO corrections 


In this paper we limit the consideration to the 
calculation of the NLO corrections to leptonic 
variables for one-photon exchange in electron- 
nucleon scattering. This approach is widely 
model independent and allows to refer to general 
non-perturbative parameterizations of the struc¬ 
ture functions which describe the hadronic ten¬ 
sor. In this way a direct unfolding of the experi¬ 
mentally measured structure functions is possible 
down to the range in Q 2 and x in which partonic 
approaches fail to provide a description of struc¬ 
ture functions. The radiative corrections calcu¬ 
lated are thus valid as well for inclusive diffractive 
ep-scattering, see e.g. [ Here the k -th order 
cross section is denoted by, 


a i = 


with C(°’°(yh, Qf) = d 2 a°/dy h dQ 2 . The O [(aX)] 
and O [(aL) 2 ] corrections were calculated in 
Ref. [ [lj[. The term (yh,Q 2 ) was derived 
in Ref. [ |] completing the 0(a) corrections. We 
re-calculated these corrections and agree with the 
previous results. 

The NLO-correction C^ 2 ’ 1 '<(y hl Q‘j) can be ob¬ 
tained representing the scattering cross section 
using mass-factorization. Although the differ¬ 
ential scattering cross section does not contain 
any mass singularity, one may decompose it in 
terms of Wilson coefficients and operator-matrix 
elements being convoluted with the Born cross 







section. In this decomposition both the operator 
matrix elements and the Wilson coefficients de¬ 
pend on the factorization scale y 2 . One writes 
the scattering cross section as, see also [ 


d 2 CT 

dyhdQ 2 


d 2 a° 

dyhdQ 2 


® Y r ei ® ® T je 

ij 


( 7 ) 


with vY(z, y 2 /m 2 ) the initial and final state op¬ 
erator matrix elements and a k i{z,Q 2 /y 2 ) the re¬ 
spective Wilson coefficients. (g> denotes a con¬ 
volution, which depends on specific rescalings of 
the chosen kinematic variables for the differential 
cross sections. Both the operator matrix elements 
and the Wilson coefficients obey the representa¬ 
tions 

r \f = 5(1 - z) + Y & m r I i f ( - rn,n) I/™”") (8) 

m>n 

a kl =S(l-z)+Y a m a^ n) L^ m - n \ (9) 

m>n 


where a = a/(2ir) and the sequences { ij} and 
{kl} in the above do always denote j(l) for the 
incoming and i(k) the outgoing particle, and L, 
L denote In (y 2 /m 2 ), In (Q 2 /y 2 ) respectively. As 
the differential cross section is /r-independent, the 
cross section is expressed by convolutions of the 
functions T ^ F ^ m,n \z) and such, that 

the /.^-dependence cancels and a structure like 
in Eq. ([|) is obtained. The present treatment in 
the OMS scheme assumes that the light fermion 
mass, m e , is kept everywhere it is giving a fi¬ 
nal answer in the scattering cross section if com¬ 
pared to the large scale Q 2 , i.e. the only terms 
being neglected are power corrections which are 
of O [(m 2 /Q 2 ) fe ] , k > 1 and therefore small. The 
last step is necessary to maintain the anticipated 
convolution structure which, in parts, is of the 
Mellin-type, as also in a massless approach. 

In the subsequent relations we make frequent 
use of the rescaling (||,||). For this purpose we 
introduce the following short-hand notation for a 
rescaling a function F(y, Q 2 ) 

Fi,f{v,Q 2 ) = F [y = Vi,FiQ 2 = Q 2 i^f) > ( 10 ) 

4 Very recently also the electron energy spectrum in muon 
decay was calculated [ E0[ using this method. 


where J, F label the respective type of rescaling. 
The NLO-corrections may be grouped into the 
following contributions : 


i LO initial and final state radiation off 

C&’ 1] (y,Q 2 ) 

ii coupling constant renormalization of 

de’%,Q 2 ) 

iii LO initial state splitting of P ie at 

C^\y,Q 2 ) 

iv LO final state splitting of P e7 at C^J\y, Q 2 ) 

v NLO initial and final state radiation off 

d°e°\y,Q 2 ) 

The function C^ 2 ’ 1 ) (y, Q 2 ) is given by 


CV’ 1 \y,Q 2 ) = Ycl 2 ' 1 \z,y h ) . ( 11 ) 


The contribution O') 2 ’ 1 ' 




= / dzP" 


6(z — 


+ dzPz. 


6{z - 


( y,Q 2 ) is 

( 12 ) 

4 )- C (M) 

z*)J F &^ - C , ( 1 > 1 )"| , 


JO J 

where C'*- 1,1 -* (y, Q 2 ) denotes the non-logarithmic 
part of the O(a) correction [ | 0 and P° e (z) is 
the fermion-fermion LO splitting function 

Pee(z) = . ( 13 ) 

1 — z 

Also the LO off-diagonal splitting functions 


P°ejz) 

P°e(z) 


Z 2 + (1 - zf 

1 + (1 -z) 2 
z 


( 14 ) 

( 15 ) 


occur in other contributions to C*- 2,1 ^. Here both 
for LO and NLO splitting functions for equal par¬ 
ticle transitions we write the contributions for 
z < 1 and account for the +-functions in explicit 
form below. 











We express the final result in terms of a{m 2 ) 
and do therfore rewrite the coupling constant by 


a(fi 2 ) 


a(ml) 




(16) 


with Po = —4/3. Due to this C^ 1,1 ' > receives the 
running coupling correction 


jl,PS,OM jIq( 0 , 0 ) 


ee,S 


■)l,NS,OM 


pl,PS,OM rFp(0,0) 
<se,T J 


e(z- z F ) j F c { p' 0) - 

( 24 ) 


C i ( i 2 ’ 1) (y,Q 2 ) = -^C( 1 > 1 )(j/,Q 2 ) . 


(17) 


The contributions Q 2 ) refer to two 


new O(a) cross sections : d 2 a ie ’^ / dydQ 2 and 
d 2 a ery ’^/dydQ 2 . The corrections are 


C§ 1 \v,Q 2 ) = / dzP^z^WC^M 2 ) 


(18) 


C^ 1 \y 1 Q 2 )= I F dzP^(z)J F (z)C^(y,Q 2 ) . 
J z o 

(19) 




The O[o 

:) sub 

Ref. 


Cpo-nM 

a 

dydQ 2 

27T 

d 2 a e T'-W 

a 

dydQ 2 

27T 


n—0,1 


( 20 ) 


= 5 : £ 

n—0,1 


( 21 ) 


Here, the functions C-j’°\y,Q 2 ) are given by 
C^\y,Q 2 ) = f dzP^(z)J I (z)C^ 0 \y,Q 2 ) 

Jz o 

( 22 ) 

C^°\y,Q 2 )= I"* dzP® e (z)J F (z)Cp’°\y, Q 2 ) . 

J z 0 

(23) 

The contribution C^" l \y,Q 2 ) reads : 

C^\y,Q 2 ) 

e( z - zl) J 1 cf ,0) - c i0 ’ 0) 


= / P 


l,iVS,OM 


ee,S 


The splitting functions in the OMS are obtained 
from the MS-splitting functions [ 0 by 


P, 


l.JVS.OM 


(*) = PefsT 5 ^ + T r - S ’ T (^) ’ ( 25 ) 


ee,S,T 


where 


^0,S,T 


(*) = -2 


1 + z 2 f 1 


(26) 


and PeS%° M (z) = PeST 5 ^)- 

the calculation are given in [ 0 - 


The details of 


4. Conclusions 

We calculated the 0(a 2 L) leptonic QED correc¬ 
tions to deep inelastic ep scattering for the case 
of mixed variables. The corrections are given in 
terms of double-differential distributions to be 
compared to the double differential Born cross 
section. The calculation was performed using the 
renormalization-group decomposition of the 2 - 
loop corrections to the differential cross section 
w.r.t. mass factorization in the OMS scheme 
for the light fermion mass. By this method an 
artificial factorization scale y 2 is introduced on 
which the physical cross section does not depend. 
Its elimination leads to a re-organization of the 
cross section which allows to assemble it in terms 
of pieces which can be calculated first individ¬ 
ually. We grouped the NLO correction into five 
terms : the LO ISR and FSR radiation correction 
of the non-logarithmic 0(a) contribution, a re¬ 
spective term due to charge renormalization, two 
new terms containing e — 7 — e initial and final 
state transitions, and the ISR and FSR OMS- 
NLO radiation correction to the Born term. 
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